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We provide a theoretical study of the quantum adiabatic evolution algorithm with different 
evolution paths proposed in [1]. The algorithm is applied to a random binary optimization 
problem (a version of the 3-Satisfiability problem) where the n-bit cost function is sym- 
metric with respect to the permutation of individual bits. The evolution paths are produced, 
using the generic control Hamiltonians H(r) that preserve the bit symmetry of the under- 
lying optimization problem. In the case where the ground state of H{ 0) coincides with 
the totally-symmetric state of an n-qubit system the algorithm dynamics is completely de- 
scribed in terms of the motion of a spin-n/2. We show that different control Hamiltonians 
can be parameterized by a set of independent parameters that are expansion coefficients of 
H(t) in a certain universal set of operators. Only one of these operators can be responsible 
for avoiding the tunnelling in the spin-n/2 system during the quantum adiabatic algorithm. 
We show that it is possible to select a coefficient for this operator that guarantees a poly- 
nomial complexity of the algorithm for all problem instances. We show that a successful 
evolution path of the algorithm always corresponds to the trajectory of a classical spin-n/2 
and provide a complete characterization of such paths. 

PACS numbers: 03.67.Lx,89.70.+c 


I. INTRODUCTION. 

Recently a novel paradigm was suggested for the design of quantum algorithms for solving 
combinatorial search and optimization problems based on quantum adiabatic evolution [2]. In the 
quantum adiabatic evolution algorithm (QAA) a quantum state is closely following a ground state 
of a specially designed slowly time- varying control Hamiltonian H(t). At the beginning of the 
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algorithm the control Hamiltonian H (0) = Hp has a simple form with a known ground state that 
is easy to prepare, and at the final moment of time it coincides with the “problem” Hamiltonian 
Hp which ground state encodes the solution of the classical optimization problem in question 
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Here E z is a cost function defined on a set of 2 n binary strings z = {zi, ... ,z n } Zk = 0, 1, each 
containing n bits. The summation in (1) is over the 2 n states |z) forming the computational basis 
of a quantum computer with n qubits. State \z k ) k of the k-th qubit is an eigenstate of the Pauli 
matrix a z with eigenvalue 1 - 2z k ± 1. If at the end of the QAA the quantum state is sufficiently 
close to the ground state of Hp then the solution to the optimization problem can be retrieved by 
the measurement. 

It has been shown recently [7] that the query complexity argument that lead to the exponential 
lower bound for the unstructured search [8] cannot be used to rule out the polynomial time solution 
of NP-complete Satisfiability problem by the quantum adiabatic evolution algorithm (QAA). 

A set of examples of the 3-Satisfiability problem has been recently constructed [4, 7] to test 
analytically the power of QAA. In these examples the cost function E z depends on a bit-string z 
with n bits, z = {z 1 ,z 2 , . . . , z n }, only via a Hamming weight of the string, w z = zi+z 2 + . . . + z n , 
so that E z = / (w z ) where the function f(w) is in general non-monotonic and defines a particular 
instance of this “Hamming Weight Problem” (HWP). In [4, 7] the original version of QAA [2] 
was applied to the HWP where the control Hamiltonian is a linear interpolation in time between 
the initial and final Hamiltonians. 

In this case, it was shown [4, 7] that the system can be trapped during the QAA in a local 
minimum of the cost function for a time that grows exponentially in the problem size n. It was 
also shown [4] that an exponential delay time in the quantum adiabatic algorithm can be inter- 
preted in terms of the quantum-mechanical tunnelling of an auxiliary large spin between the two 
intermediate states. 

The above example has a significance greater than just being a particular simplified case of a 
binary optimization problem with symmetrized cost. Indeed, one can argue that it shows a generic 
mechanism for setting “locality traps” in the 3-Satisfiability problem [10], But most importantly, 
this example demonstrates that exponential complexity of QAA can result from a collective phe- 
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nomenon in which transitions between the configurations with low-lying energies can only occur 
by simultaneous flipping of large clusters containing order-n bits. In spin glasses, there is typically 
an exponential number of such configurations, the so-called local ground states. A similar picture 
may be applicable to random Satisfiability problems [15]. In some cases, these transitions can be 
understood and described in terms of macroscopic quantum tunnelling. A tunnelling of magne- 
tization was observed in large-spin molecular nanomagnets [11] and in disordered ferromagnets 
[ 12 ]. 

The paper [1] suggests that large tunnelling barriers can be avoided in QAA by using multiple 
runs of QAA with realizations of the control Hamiltonians H(t) sampled from a random ensem- 
ble. This ensemble is chosen in a sufficiently simple and general form that does not depend on the 
specific instance of the optimization problem. Different Hamiltonians H{r) correspond to differ- 
ent paths of the unitary evolution that begin and end in the same initial and final states (modulus 
phase factors). The complexity of QAA with different paths for the HWP was tested numerically 
in [1] using an ensemble of random 8x8 matrices. The results indicate that the HWP may be 
solved in polynomial time with finite probability. 

In case when the random paths H{r) preserve the bit-permutation symmetry of the problem 
it is natural to describe the random ensemble of H(t) in terms of the dynamics of a spin-n/2 
system. This approach allows for a general theoretical analysis of the algorithm. In the present 
paper, we perform this analysis for the random version of HWP (over-constrained 3-Satisfiability 
problem) by mapping the dynamics of QAA onto the motion of a quantum particle in a ID effective 
potential. This allows us to compute the statistical weight of the successful evolution paths in the 
ensemble and provide a complete characterization of such paths. 

II. QUANTUM ADIABATIC EVOLUTION ALGORITHM WITH DIFFERENT PATHS 

In a QAA with different paths [1], one specifies the time-dependent control Hamiltonian 

m = h( T ) 


H(r) = (1 - t)H b +t(1 -t)H e + tH f> (3) 

r = f € ( 0 , 1 ). 


where the control parameter r plays the role of dimensionless time. This Hamiltonian guides the 
quantum evolution of the state vector \ip(t)) according to the Schrodinger equation ihd\'ilj(t))dt = 
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i?(r)| , 0(t)) from t = 0 to t = T, the run time of the algorithm. Hp is the “problem” Hamiltonian 
given in (1). H B and Hp are ‘driver” Hamiltonians designed to cause the transitions between the 
eigenstates of Hp. 

An initial state of the system [0(0)} is prepared as a ground state of the initial Hamiltonian 
77(0) = Hp- It is typically constructed assuming no knowledge of the solution of the classical 
optimization problem and related ground state of Hp. In the simplest case 

n 

= h/’(0)) = 2-“/ 2 X>), (4) 

3 = 1 z 

where a J x is a Pauli matrix for y-th qubit and C > 0 is some scaling constant. The ground state of 
Hp has equal projections on any of the 2 n basis states \z) (2). 

Consider instantaneous eigenstates |<fo.(r)) of H(t) with corresponding eigenvalues E k (r) ar- 
ranged in non-decreasing order at any value of r 6 (0, 1) 

H{r)\ ( p k {T)) = \ k {r)\ ( t>k{r)), k = 0, 1, . . . , 2 n - 1. (5) 

Provided the value of T is large enough and there is a finite gap for all t 6 (0, T) between the 
ground and exited state energies, AA(r) = Ai(r) — Ao(r) > 0, quantum evolution is adiabatic 
and the state of the system stays close to an instantaneous ground state, |0 o (i/T)) (up to a 
phase factor). Because H(t ) — H P the final state [0(T)) is close to the ground state | 0 o (r = 1)) 
of the problem Hamiltonian. Therefore a measurement performed on the quantum computer at 
t — T will find one of the solutions of combinatorial optimization problem with large probability. 
Quantum transition away from the adiabatic ground state occurs most likely in the vicinity of the 
point T rti T c where the energy gap AA(r) reaches its minimum (avoided-crossing region). The 
probability of the transition is small provided that [17] 

T>hH max AA"? nJ (6) 

where 

Hxnax = m . ax . I (01 ( r ) I < ~J~ 1 00 ( r )> ! > 
t€(o,i) dr 

AA min = min [Ai(r) - A 0 (r)] , (7) 

TS(0,1) 

The r.h.s. in Eq. (6) gives an upper bound estimate for the required runtime of the algorithm and 
the task is to find its asymptotic behavior in the limit of large n » 1. The numerator in (6) is of 
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the order of the largest eigenvalue of dH/dr = H P — H B + (1 — 2t)H e , which typically scales 
polynomially with n. However, AE min can scale down exponentially with n and in such cases the 
required runtime of the quantum adiabatic algorithm to find a solution grows exponentially fast 
with the size of the input. 

One should note that the second term in the r.h.s. of (3) is zero at r = 0 and r = 1. Therefore, 
by using different driver Hamiltonians H E one can design a family of (possibly random) adiabatic 
evolution paths that start at r = 0 in the same generically chosen initial state and arrive at the 
ground state of Hp at r = 1. In general, different paths will correspond to different mi nimum 
gaps g m i n and one can introduce the distribution of minimum gaps. This distribution can be used 
to compute the fraction of the adiabatic evolution paths / that arrive at the ground state of H P 
within polynomial time, 

T < cn~ a , a > 0, c=0( 1). (8) 

For a successfully designed family of paths the fraction / is bounded from below by a polynomial 
in 1/n which leads to the average polynomial complexity of QAA. 

ID. BINARY OPTIMIZATION PROBLEM WITH SYMMETRIC COST FUNCTION 

Consider a binary optimization problem defined on a set of n-bit strings z with the cost function 
E z in the following form: 

n 

E 2 = f (w z ) , w z = Y z r ( 9 ) 

i = i 

This cost is symmetric with respect to the permutation of bits, it depends on a string z only through 
the number of unit bits in the string w z (the Hamming weight). In this paper we consider the cost 
function (9) in the following form which is generalization of the cost introduced in [4, 7, 9] 

E z = c ( z h+ z i2+ z i 3 )> ( 10 ) 

il <^2 <^3 

c(m) = PoSmfi + Pi5m : i + P2Sm,2 + P35 m ,3- 

Here the sum is over all possible 3-bit subsets of the n-bit string z. A subset z ix + z i2 + z i3 
contributes to the total cost a weight factor p k where A: is a number of units bits in the subset. A set 
of weights {pk} defines an instance of this generalized Hamming Weight Problem (HWP). One 
can formulate a random version of HWP, e.g., by drawing numbers {p k } independently from a 
uniform distribution defined over a certain range. 


6 


In the limit of large n » 1 the cost function (10) takes the following form: 

3 

E z = l 3 Gp( 1-y), Gp(q) = '£&<?, (11) 

fc =0 

here l = nj 2 and we only keep the terms of the leading order in n. The coefficients (3 k in (11) are 
linear combinations of p k 

Pk = y [pi + (— 1 )^ 2 ] + g [po + ( 1 )^ 3 ] • (12) 

here £* = 1 for k = 0, 1 and £ k = — 1 for k = 2, 3. 

The function Gp(g) in (11) is a third degree polynomial in q, and the form of the function 
depends on the coefficients /3 k (Pk)- H I s eas y t0 show that there is a finite size region in the 
parameter space {p k } where Gp(q ) is a non-monotonic function of q that has global and local 
minima on the interval q € (— 1,1). Those minima are separated by a finite barrier with width 5q = 
0(1). The barrier separates strings that have close values of the cost E z but are at large Hamming 
distance from each other: they have 0(n) distinct bits. This property can lead to exponentially 
small minimum gaps in QAA due to the onset of low-amplitude quantum tunnelling [4]. 


IV. CONSTRUCTION OF THE CONTROL HAMILTONIAN 

A. Representation in terms of operator components of the total spin 

It is natural to consider the control Hamiltonians (3) for solving the HWPs that are symmet- 
ric with respect to permutation of individual bits (2). In what follows, we use the normalized 
components of the total spin operator S for the system of n individual spins- 1 



Here Sj are the projections of the total spin operator on the j-fh axis (j = x, y, z) and <f] are Pauli 
matrices for the i-th spin. For the sake of bookkeeping, in (13) and also throughout the paper we 
use “hats” for the spin operators, such as Sj, hj, and some others, in order to distinguish them 
from their correspon din g eigenvalues (Sj and rij, respectively, in the above example). 

To obtain the problem Hamiltonian (1) we make use of the obvious connection between the 
values of the H amming weight function w z of an n-bit string z and corresponding eigenvalues n z 
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FIG. 1: Plots of the cost function (11) G p vs q for different choices of the weights {pk}- Curve 1 corresponds 
to po = 0, pi = 3, p 2 =- 1, P 3 = 1, and the cost function Gp{q) has a global minimum at q = 1, 
corresponding to the string z with the Hamming weight zero, z\ = — . . . — z n = 0. It also has a local 

minimum at q = — 1 corresponding to the bit string with Hamming weight n, z\ = zi = . . . = z n = 1. The 
curve 1 yields the particular form of the cost function Gp(q) considered in [4], [1], Curve 2 corresponds to 
po = 0.5, pi = 2.5, P 2 = —2, pz = 0.3, it has a global minimum at q — q* inside of the interval (—1, 1). 
This minimum corresponds to approximately ( n ^.) bit strings z that all have the same Hamming weight 
w z = w* = n(l — q*)/ 2. 

of the spin projection operator n z 

7 n 

n 2 |z) - n 2 |z), n z = 1 — (14) 
Then from Eqs. (2), (11) and (14) we obtain 

H P (n z ) = l z Gp{n z ). (15) 


We chose the driver H B in a bit-symmetric form that coincides with (4) (up to a constant term) 

H b = l 3 G B (n x ), G b (x) = 2(1 - x). (16) 
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B. Bit-symmetric drivers H E 

It was proposed in [1] that H E can be constructed using some generic ensemble of random 
matrices. The bit-symmetric random drivers for the cost functions of the type (10) can be con- 
structed as follows [1]. One generates an 8 x 8 random Hermitian matrix A with zero diagonal 
elements and non-diagonal elements that are independent random numbers identically distributed 
in a certain interval. Matrix elements of A ZiiZjiZk can be enumerated by all possible configurations 
of a 3-bit string {z u Zj, z^}. Then H E takes the form 

E ww, (i7) 

i<j<kze{ 0,l} n 

Here \z) are computational basis states Eq. (1) corresponding to bit-strings z = {zi, . . . , z n }. Each 
randomly selected A generates H E and therefore a random path modification of the QAA. 

From the above discussion, it follows that the matrix of the operator H E (17) is symmetric 
with respect to the bit permutations and therefore it commutes with the operator of a total spin 
S 2 of a system of n spins This means that H E acts independently in each of the sub-spaces 
corresponding to certain values of the total spin 0 < / < f [22], It follows from (15) and (16) that 
the same is true for the total control Hamiltonian (3) 

[ff(r),S 2 ] =0, r € (0, 1). (18) 

Since in our case the initial state (4) is a totally symmetric combination of all states and therefore 
corresponds to the maximal spin l = |, our system always stays in this sub-space during the 
algorithm. Therefore in the analysis of the complexity of QAA one can reduce the 2 n x 2 n matrix 
of H (r) to the (2n+ 1) x (2n+ 1) matrix that only involves the states with different spin projections 
of the maximum total spin l = |. Binary strings corresponding to the quantum states from this 
subspace are distinguished from each other by their Hamming weight only. 

In Appendix A, we show that in the case of real- valued symmetric matrices A and in the large- 
spin limit, the bit- symmetric driver He (17) can be presented as a linear combination of 6 operators 
expressed in terms of the large spin operator components n x , h z acting in the subspace with l — 
Using this fact, and also Eqs. (15) and (16) one can write a bit-symmetric control Hamiltonian (3) 
in the following form 


jzH(r) = G(t, n x , n z ) = (1 - r) G g (h x ) + r(l - r) G B (n x , h z ) + tG p (n z ) , (19) 
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G E (n x , n z ) = 71 n x + 72 n 2 + 73 h\ + 74 n x n z + 7 5 n x n 2 + 7 &n 2 x n z , (20) 

where {7*} ( k — 1 , . . . , 6) are independent real coefficients given in Eq. (A8). As we show in 
Appendix A, any random realization of the matrix A can be mapped onto combinations of drivers 
(20) by the appropriate choice of the coefficients 7*. This fact allows us to analyze the minimum 
gap in QAA with different paths (3) using the WKB analysis of the dynamics of a spin-| in the 
large spin limit (n » 1). 

V. ADIABATIC EVOLUTION OF A LARGE SPIN 
A. WKB approximation for the large spin 

Our analysis in this section is a particular application of the WKB -type approach commonly 
used for the description of quantum spin tunnelling in magnetics [18], [19], [20], [12], This 
approach is applicable for the large spins (l » 1), which is the case of interest for us. 

We choose a as a quantization axis and following the standard procedure to obtain the effective 
quasi-classical Hamiltonian in polar coordinates { 6 , 7?} with 9 G [0, 7r] and p G [0, 2-7r]. We make 
use of the Villain transformation [23] 

n x = \fl + e - n z (n z + e) cos {(p ) , e = J> 

where azimuthal angle operator (p satisfies the commutation relation 

[<P, fiz] = ( 22 ) 

In a change of notation we introduce a coordinate q and canonically-conjugate momentum p (cf. 
[23]) 

ci 

q = n z , p= -ie— = - <p , (23) 

dq 

(—1 < q < 1). Expanding (21) in the large spin limit e <§C 1, we obtain 

n x = (1 - <7 2 ) 1/2 cosp + + °( e2 )- C 2 4 ) 

Finally, we write the scaled Hamiltonian of the system (19) in terms of the new variables 


G(n x ,n z ,T ) = 'H(q,p,r), 
ft(?,P,r) = K(q,p,r) + AH(g,p,r) 


( 25 ) 
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where 


H(q,p,r) = G (yi - g 2 cosp,g,r) , 

[q,p] = 


and AH is a small correction 


ah (^> = ^£ +0 < £2) 
(here dG/dn x has the same arguments as G in (26)). 

The stationary Schrodinger equation (5) in the new basis 

n(q,p,r)^ k (q;T) = A fc (r)'I' fc (q; r), 


(26) 


(27) 


(28) 


(29) 


can be solved in the WKB approximation with the small parameter e < 1 playing the role of a 
Plank constant. Then the wave function \&(#) takes the form 

'iA{q)' 


V(q) = B{q) exp 


(30) 


where in the leading order in e the action function A(q) satisfies the following Hamilton- Jacobi 
equation 

nfq,- ^,r)=A. (31) 


dq 

This equation describes a ID auxiliary mechanical system with coordinate q, momentum p, energy 
A, and Hamiltonian function H(g, p, r). Classical orbits satisfy the Hamiltonian equations 


q(t) = H P (q(t),p(t),r), p(t) = -E q (q(t),p(t),r), 


(32) 


where H q and H v stand for the partial derivatives of H with respect to q and p, respectively. Sta- 
tionary points of the dynamics (q*,p*) correspond to the elliptic and saddle points of the Hamilto- 
nian function 

H 9 0*,p,,r) = 0, H(?*,j>„t) = 0, (33) 

Elliptic points are minima (or maxima) of H (q, p, r) on the (q. p) plane. They satisfy the condition 

All = Hpp(q*,p*, T)H qq (q*,p*, r) - H qp (q*,p*,r) > 0. (34) 

where H PP is understood as a second derivative of H with respect to p, etc. Saddle points corre- 
spond to Ql < 0 in (34). 
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In the limit e <C 1 the adiabatic ground state 'I ' 0 (g; r) (29) is localized in the small vicinity of 
the fixed points (g*,p*) corresponding to the global minimum ofH(g, p,r) atagiven value of r. To 
logarithmic accuracy the WKB-asymptotic (30) of the ground-state wave function is determined 
by the mechanical action for the imaginary-time instanton trajectory (q(t),p(t)) emanating from 
the fixed point (q t , p* ) 


^0 (?;r) « # o(q',q*,P*,T ) oc exp 


% 

U 


o 


dtq(t)p(t ) 


q(-ioo) = g*. p(-ioo) = p», g(0) = g. 


(35) 


Integration in (35) is along the imaginary axis (— ioo, 0). The instanton trajectory obeys Eq. (32) 
with the boundary conditions given above and f 6 (—zoo, 0) corresponding to the line of integra- 
tion in (35). The choice of the final instant, t = 0, is arbitrary since the instanton trajectory is 
degenerate with respect to a shift of the time axis. 

We note that the WKB asymptotic (35) decays exponentially fast as the coordinate q in (35) 
moves away from its value at the global minimum g* into the classically inaccessible region. This 
corresponds to the growth of the imaginary part of the action in (35), similar to the conventional 
quantum tunnelling in the potential. In the vicinity of (g*,p*) the ground-state wave function 
&o(q) takes the form similar to that of harmonic oscillator: 


Mq) 


m* = 


c x exp 
1 


- (p t 5q — 5q 2 ) - 

e V q 2 H pv q J 2e 


ITT ( "TT ) (36) 

I Hp P {q*,P*>T)\ 

here cu, > 0 is defined in (34). Similarly, the energy spectrum in that region corresponds to the 
classical elliptic orbits with oscillation frequency 0* 


A*. - H (q*,p*,r) ~ e£l t 


k+ l 


k = 0,1,.... 


(37) 


We note that the frequency fl* depends on r and determines the time- varying instantaneous gap 
between the ground and first exited states, A A = ef2»(r). 

VI. LOCAL AND GLOBAL BIFURCATIONS DURING THE QAA 


It can be seen from Eq. (26) that the global minimum of H(g, p, r) will correspond to p* = ±kir 
(k = 0, ±1, . . .) provided that the following condition holds for all n x : 

_ <9G(n x ,g*(r),r) 


OOt 


dn r 


0, 


(38) 
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where the positive and negative signs of u x correspond to even and odd values of k, respectively. 
The value of q* in (38) corresponds to the global minimum of the effective potential U(q, r) 


U( q ) = G(v / 1^7,g,r), U(q)-U[q.)> 0. (39) 


Under the above condition the Hamiltonian function of the system near the global minimum (g*, 0) 
exactly corresponds to that of the harmonic oscillator with effective frequency (34) and mass 
m* (36) 


H (q,p,T) = 


2m* (r) 


V + 


m*Ql(r)(q - q*{r)f 


- = -/r r 7i 

m* 

m*Ql = U"(q*). 


dG(y/l-qlq*,r 


dn r 


(40) 


In the WO picture the ground state of the system correspond to the particle performing zero-level 
oscillations near the bottom of the slowly varying potential U(q,r). There are two types of the 
bifurcations that can destroy the above adiabatic picture: 


1. Local bifurcation 


Assume that at some instant of time r = r 0 the effective mass m* (r) goes to infinity. In the 
vicinity of this point the Hamiltonian function (26) can be approximated as follows: 


H (q,P, r) = ~p 4 - ^sp 2 + ^ Sq 2 + d 0 sSq + 0(s 5/2 ), 
5q = q- q*(r Q ), s = t - r 0 , 


(41) 


where 


o 0 = 


Co = 


d 2 G n 2 ( u , d 2 G r, 

(i - q*( T o)), b o = - q*( T o), 


dn 2 x 

d 2 U 


dn x dr 


do = 


d 2 U 


(42) 


dq 2 ’ 8qdr 

in the above equations all functions are evaluated at the point (q*(ro),p = 0). Equation (41) 
corresponds to A 3 bifurcation point [24], It can be seen from (41) that for r > r 0 the single global 
minimum of H(q, p, r) splits into the two minima with nonzero momenta 


±, \ ^ , /'6&o(t-t o)V /2 ^ n 

PI (t) ± , r - r 0 > 0. 

&o J 


( 43 ) 
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Due to the symmetry H(q,p, r) = H(g, -p, r) the two global minima with nonzero p* will stay 
symmetric with respect to the #-axis at later times. 

It follows from (34), (40)) that the linear oscillation frequency vanishes at the bifurcation point, 
Q*(r 0 ) = 0, however the energy gap AA (r 0 ) ^ 0. By solving the Schrodinger equation (5) at this 
point in the representation of the momentum p one can find the eigenfunctions 'IQ (p, r 0 ) and eigen- 
values Afc(ro) corresponding to a ID quantum system moving in a quartic potential (cf. Eq. (41)). 
This analysis yields an estimate for the value of the gap, and the characteristic localization range 
Sp for y k (p,T 0 ) 

AA ~ e 4//3 , Sp ~ (44) 

The size of the energy barrier in momentum p separating the two global minima in (41) grows 
with time for r - r 0 and this leads to a rapid decrease of the energy gap. Sufficiently far from 
the bifurcation point, r - r 0 » e 2/3 , each of the global minima (q*(r),p^(r)) gives rise to its 
own WKB asymptotic (35) localized at the minimum. The ground state and the first exited state 
correspond to their symmetric and anti- symmetric combinations, respectively 

Wk(q) = -4= (y(q;q*,pt,T) + (- 1 )^( 9 ; q*,v~,r )) , 

k = 0, 1. (45) 

For r ~ t q e 2 / 3 the tunnelling splitting of energy levels for the symmetric and antisymmetric 
states determines the value of the gap AA (r) and decreases exponentially fast with r — tq. Away 
from the bifurcation region, r - r 0 = 0(1), the gap scales down exponentially with n (note that 
e = 2 )n). 

As a result of the local bifurcation, the purely adiabatic evolution in QAA collapses. The 
amplitude of staying in the adiabatic ground state for r < r 0 is nearly equally split between the 
states (45) with the two lowest eigenvalues. In general, this may reduce the probability of finding 
a system in a ground state at r = 1 by a factor of 2. We note that the control Hamiltonian (19) 
is at most a cubic polynomial in n x , n z , and therefore the number of local bifurcation events 
during QAA is of the order of one. In the worst case they will cause the reduction of the success 
probability in QAA by a constant factor. 

For a given instance of the cost function (11) defined by the coefficients 0 k (or p k ) the onset of 
local bifurcations (41) depends on the choice of the driver Hamiltonian H E (20). 

There are a number of ways to select coefficients y^’s in the driver Hamiltonian (20) to avoid 
local bifurcations during QAA in a broad range of values of the coefficients p k . For example, to 
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completely suppress local bifurcations (41) one can 'keep in (20) only terms linear in h x and set 

72 = 73 = 76 = 0. (46) 


A. Global bifurcation 


The Hamiltonian function H = H(g,p , r) defines a 3D surface over a 2D plane (g,p) and the 
shape of this surface varies with time r. We consider global bifurcations of this surface where 
the energies of its two minima cross each other at some instant of time t — tq while the distance 
between the minima on the (q,p) plane remains finite at the crossing point. For r > r 0 the 
minima exchange their roles: global minimum becomes local and vise versa. Before and after the 
intersection in the energy space the two minima are uniquely identified with the ground and first 
exited states of the system’s Hamiltonian (26). The corresponding wave functions # 0,1 (?) are well 
approximated by their asymptotic expressions(35),(36). 

The small vicinity of the global bifurcation point can be described within the standard 2-level 
avoiding-crossing picture. There # 0,1 (?) are given by symmetric and antisymmetric superpositions 
of the WKB -asymptotic corresponding to intersecting minima. The value of the gap changes with 
t im e as y / c 2 (r — To) 2 + AA^ in where c is some constant and the minimum gap is determined by 
the overlap of the WKB asymptotic. To logarithmic accuracy it is given by the imaginary part of 
the mechanical action (35) along the instanton trajectory connecting the two minima 


6 log AA m i n 
lim q{r) = 

T—>±2 CO 


Im 


dt q(t)p{t) 


9*’ 2 > lim p(t)=pI’ 2 

T—+±lOO 


(47) 

(48) 


Here q,,p* are coordinates of the two minima; H(ql,pl,T 0 ) = H{q 2 ,pl, r 0 ), and the instanton 
trajectory obeys the Eqs. (32). The analytical expression for the minimum gap was studied in [4], 
[9] for the case H E — 0, using a simplified version of the Hamming Weight problem (10). Below 
we identify certain geometrical properties of the global bifurcations in the case He = 0 that will 
be used later in the selection of the drivers H E for the successful QAA. 
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FIG. 2: The global bifurcation mechanism: the effective potential profiles U (q, r) vs q for r < To, r = tq 
and t > tq are represented by the curves 1,2, and 3, respectively. 


1. The case He — 0 


In the case = 0 (j = 1, ...6), the Hamiltonian has a minimum at p* (r) = irk and the value of 
g*(r) corresponds to the global minimum of the effective potential U(q, r) (39). We use Eq. (19) 
and also the Condition U'(q*) — 0 to obtain the following equation for q*(r) 

d Q*( T ) G ffejr)) 

dr (1 -T)U»(q.(T),Ty 1 ' 

This equation holds until the global bifurcation point at r = r 0 where q*(r) changes discontinu- 
ously in time (see Fig. 2). At the minimum of the potential U"(q ,) > 0 and therefore the direction 
of the motion of q*(r) entirely depends on the direction of the “force”, — G" P (q*). At r = 0 the 
potential U(q, 0) has a unique minim um at the point q = q*(0) = 0. It is clear that with this 
initial condition equation (49) can lead to a “wrong” minimum of Gp{q) that lies above the global 
minimum, and such cases will give rise to a global bifurcation. This effect is illustrated in Fig.l 
where the two different cost functions correspond to the same direction of motion for q*{j). The 
value of q M (t) may either smoothly approach the global minimum of the cost (curve 2), or move 
toward a ’’wrong” local minim u m (curve 1), leading to the global bifurcation and exponentially 
small gap in QAA. Adding He to the control Hamiltonian can invert the direction of motion of 
q*(r) toward the global minimum of Gp(q). This can be seen from the fact the Ea.(49) in presence 
of H e possesses the additional term 

(1 — r) dG E (y/l ~ 

U"(q*, r) dq * 


( 50 ) 



16 


(here we drop for sake of brevity the argument r in q*(r)). Clearly, the successful G E should not 
possess reflection symmetry with respect to n z . Therefore we should only select the terms in (20) 
that contain odd powers of n z . Taking into account (46) we arrive at the following form of the 
driver Hamiltonian 

G E (n x ,n 2 ) = 7 4 n 2 n z . (51) 

This driver can remove the potential barrier between the two competing global minima of U (q, r) 
by shifting the original minimum at r = 0 towards the true global minimum of the cost function 
Gp(q) (cf. Fig. 1). In the classical picture (26) the driver (51) corresponds to an external field 
parallel to a-axis which can destroy the tunnelling barrier along this direction. The mechanism of 
such tunnelling avoidance is similar to the one considered in [1], where the external field gener- 
ated by the driver (51) compensates the effective field due to the linear term proportional to the 
coefficient 0i in the problem Hamiltonian (11). 


B. Bifurcation transition to the tunnelling regime 


In general, one can expect that a complete suppression of the tunnelling barrier at all values 
of r requires a certain magnitude (and sign) of the coefficient 74 depending on the choice of the 
coefficients @ k in the cost function Gp{q). 

The transition to the tunnelling regime can be described as an A 3 bifurcation point, illustrated in 
Fig.2. The effective potential U changes parametrically with r, 74 and {j3 k }. Near the bifurcation 
point (r c , 7 4c) q c ) ,the potential has the form U = a 5q 4 + b SjSq 2 + cSq5r where Sr, 6 q, 5-y are 
deviations from the bifurcation point in r, q and 74, respectively. The corresponding conditions for 
the A 3 bifurcation point are: 


dU 

dq 


d 2 U d 3 U 


= 0 . 


dq 2 dq 3 

Taking into account (39) and (19), (20), (51), the above equation yields 

_ n , 2 t c 2 ( 1 - t c )( A +/? 3 ) 2 

W-fh) = 1 + t to To — TTI2 ’ 

3 [2 - r c (2 - #>) ] 

_ (1 ~ r c) (3/?i + Pz) 4- t c (3 2 P 3 

74c f J- T cj — n o /-> _ \ > 

T c p2 2(1 T c ) 


(52) 


(53) 


These equations should be solved for 7 4c and r c for the given set of the coefficients fit-- The 
bifurcation is avoided when 


M > |74c|- 


(54) 
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For example, in the particular case of the HWP (10) considered in [4], [1], we have 

Pi = 1 / 2 , /% = -3/2, & = -7/6, (55) 

r c ~ 0.44, 74c ~ —0.95. 

In this case the example of the driver Hamiltonian H E that allows to avoidance of tunnelling in 
QAA was given in [1] where the value of 74 = —8 was used. According to (55) this value is way 
below the critical value 74 c . 


1. Numerical Simulations of the bifurcation boundary 

We performed numerical simulations with the effective potential (39) checking for the onset 
of tunnelling for all p k 6 [0; 3], k = 1, ...4. The numeric simulations confirm that the situation 
discussed above is typical for the general HWP, implying that (51) is the only driver term that can 
be fundamentally responsible for the tunnelling avoidance in a general case, if the coefficient 74 is 
defined appropriately. In particular, one of the two drivers (51) with 

74 > 7c = 4.9 or 74 < -j c = -4.9 

Pk € [0; 3] k = 0,1,2, 3, (56) 

always suppresses tunnelling in the QAA. 

16 

12 

£ 8 
4 
0 



FIG. 3: The critical value j c vs domain size L. 
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We also solve the Eqs. (53) numerically for coefficients p k taking values on a dense grid of points 
in the cube p k 6 [0; L\ (k= 0, 1,2,3). For each size of the cube L we select the point with largest 
value of 74 c denoted below as j c = 7 C ( L ). The results are presented in Fig.3. The critical value 
7 c is monotonically increasing in L, and the dependence is close to linear for sufficiently large 
L, but it is non-linear in the range 0 < L < 3. It can be inferred from Eq. (53) that a nonlinear 
dependence of 7 C on the scale L is due to the fact that the critical time r c also depends on L. 

The linear dependence of y c ( L ) for large L has a simple intuition. According to (11) and (12), 
the magnitude of Gp(q) is proportional to L. According to Eq.(49), the maximal magnitude of 
the coefficient presents a ’’force” that can possibly move a system into the local minimum at 
small r . From ( 12 ), we conclude that |Al majc = max p fc e[-z,zj Pi = 2T. In the limit of large L, the 
role of the driver G B in (53) becomes unimportant. Therefore, the only competing terms are the 
driver Ge and the problem Hamiltonian Gp. The term (50) generated by Hp compensates (50) 
the ’’force” fa when |y 4 | > \Pi\ m&x , and therefore in this limit we have 

7 c (L) = max y 4c « \Pi\ max = 2 L. (57) 

p k e[-L,L] 

One should note that among the effective potentials generated by choosing different {p k }, there 
are two subsets that can be mapped onto each other by means of the mirror reflection about the 
q-axis, U(q,r ) — ► -U(q,t). We note that the same driver He can not simultaneously suppress 
tunnelling barriers in each of the two mutually symmetric potentials: if the tunnelling barriers are 
not suppressed with j C) they will be suppressed with — 7 C , and vice versa. This gives a simple 
intuition for the tunnelling barrier suppression boundary (54). 

Finally we conclude, that it is possible to indicate the range of value of ( 74 1 such that the driver 
Hamiltonian H E = l 3 74 n x n z will play the role of a universal driver that guarantees polynomial 
performance of the QAA for all instances of the generalized Hamming weight problem (10) provi- 
sory to the mirror-reflection symmetry in the possible choice of the cost functions and the common 
normalization factor L. 

VII. PROBABILITY OF SUCCESS OF THE QAA WITH RANDOM PATHS 

Using the analysis from the previous section one can estimate the probability of success for 
the QAA with random paths proposed in [1], In that algorithm, the ensemble of random drivers 
H e was generated using random 3x3 matrices A Zi z . Zk (17). It is shown in Appendix A that 
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for the bit-symmetric optimization problem (10) the above ensemble is identical to the ensemble 
of independent uniformly-distributed random coefficients 7^ (k=l — 6) that appears in the large- 
spin representation of the driver He (20). Then for any instance of optimization problem in (10) 
defined by the set of the coefficients {; p k } one should compute the fraction / of the domain of the 
coefficients {7*.} where the following conditions are satisfied: 

(i) . Condition for the nonzero effective mass (38). 

(ii) . The condition (54) for the complete avoidance of the tunnelling barriers in combination with 

Eq.(53) for the bifurcation boundary. 

Here we compute the fraction / for the particular instance of the optimization problem (10) con- 
sidered in [1, 4, 9]. In this case Eq. (38) takes the form 

■7—^ = —2 + T [71 + 27271* + 373^ + 7 4 q 

1 — T 

+I 5 Q 2 + 275 71* q] 7^ 0, (58) 

where n x = ±y/l — q* and q* provides global minimum of U (q, r) (39). The effective mass is 
non-zero if u x 7^ 0, and (58) yields an estimate on the range of {7^} as 

I72I + l 7 e| <1 + 1/2 (I71I + 3 |7 3 | + |7 4 | + |7 5 |) . ( 59 ) 

Following [1] we assume that the non-diagonal matrix elements A ZiZ . Zk are distributed in the 
interval [—3, 3]. Making use of (A8), we obtain 

I72I + |ts| < 16 , I71I + 3 (73 1 + I74I + I75I < 50 . ( 60 ) 

Therefore, the probability that inequality (59) is satisfied is estimated as 1-15 2 / (50 x 16) ~ 
0.71875. On the other hand, the values of 74 in (A8) belong to the range, —12 < 74 < 12. Using 
the value of 7 c « -0.95 given in (55) we estimate the probability of 74 < — 7 C to be approximately 
equal to 0.46. Making an approximation that the cases when the effective mass is non-zero are 
statistically independent from the cases when 74 < — j c , we obtain the total probability of success 
as P tot se0.46 x 0.71875= 0.334 « 1/3, which is in qualitative agreement with the numerical 
results of [1]. This estimate can be generalized to the case when the matrix elements A ZiZ . Zk are 
distributed in the interval [—L,L\ for sufficiently large L > 3. In this case, the probability that 
(59) is satisfied remains the same, « 0.71875, while the probability that 74 < — 7 C is estimated as 


20 


(4 L — | 7 C |) /8L. With the assumption of statistical independence, the total probability of success 
is P tot « 0.718 x (4L — 0.95) /8L, and in the limit of large L » 1 we have P tot ps 0.359 which 
exceeds slightly the value for L = 3. 


VIII. POLYNOMIAL QAA AND CLASSICAL DYNAMICS OF LARGE SPIN 


In absence of tunnelling, the dynamics of the large spin can be characterized by classical equa- 
tions of motion for the spin projections treated as c-numbers in the form [18] 



(61) 


with 


dH _ f dH dH dH ] 

" ~<95 " W dS y ' dS 2 J- 

In coordinate form and in terms of the dimensionless spin projections, this yields 


(62) 


dt 

dn y 

dt 

dn z 

dt 


-u z n y , 

Z^x 

j 


(63) 


where we took into account that since H does not contain the S y component, To = = 

{(jj X} 0, co~}. In the case when the ’’effective magnetic field” uj does not explicitly depend on 
time, the system (61), (62) has two independent integrals of motion 


5 2 = 5* + 5j + 5 2 , 


(64) 


J = - TJ,S = 
to 


S X 4" Z 5r 


A 2 + 

The first (64) reflects the conservation of total spin and also holds for an arbitrary time-dependent 
field do = |=j), whereas the second integral corresponds to the adiabatic invariant of the system 
(61), (62). Since in our case To — is parametrically time-dependent, the adiabatic invariant 
is conserved approximately for sufficiently slow parametric evolution. Note that the adiabatic 
solutions always play the role of ’’envelope solutions”. This means that on average, the spin 
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closely follows the adiabatic solution, but there are fast oscillatory-type motions superimposed 
on the slow adiabatic evolution. Basically, the adiabatic approximation in the classical case is 
applicable when the ’’slow” motion is much slower than the fast oscillatory motion. This exactly 
corresponds to the adiabatic evolution of the spin system in the quantum case [ 22 ], 

Making use of (64) and taking into account that at the instant r = 0, the total spin was parallel 
to the z-axis, we obtain J II = = i ; or 

v Zj 

~n = — , (65) 

U! 

implying that the total spin is always parallel the effective magnetic field To. Therefore, the adia- 
batic evolution of the large spin can be simply described as the situation when the spin follows the 
effective field (on average). 

We note that at this level, there is a direct correspondence between the adiabatic classical so- 
lution and the quasiclassical wave functions of the large spin parallel to ~n. From (65), it follows 
that this direction can be identified with the effective magnetic field To = This justifies the 
’’variational” approach introduced in [1, 4], identifying the variational wave functions with the 
adiabatic ground states along the evolution paths when the total spin is parallel to fl . Therefore, 
one can observe that in the absence of tunnelling, the general HWP is solved essentially by the 
classical paths of the QAA. 


IX. CONCLUSIONS 

We apply the quantum adiabatic evolution algorithms with different paths [1] to the general- 
ized Hamming Weight Problem that corresponds to the specific case of the random Satisfiability 
problem defined in (10). We show that any random evolution path produced by this algorithm for 
the HWP can be obtained by using 6 specific deterministic basis operators with random weights 
and therefore is parameterized by 6 independent random numbers. Therefore, the approach to 
QAA with different paths can still be reduced to the large spin dynamics for the HWP. We show 
that only one of these ’’generators” can be a ’’universal” driver fundamentally responsible for tun- 
nelling suppression for arbitrary HWP and therefore the problem of constructing such a universal 
driver reduces to the definition of its weight 74 . Due to the possible reflection symmetry of the cost 
function, any particular case of the general HWP can be solved with one of the two values of the 
weight with | 7 4 | > 7 4c , that is by applying one of the two universal path modifications. We analyze 
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the nature of the wave functions along the successful paths and show that it is quasiclassical and 
corresponds to the dynamics of a large classical spin. Therefore, we show that the general HWP is 
solved by completely classical paths of the QAA and present a complete characterization of these 
paths. 

We analyzed in details the types of bifurcations of the effective Hamiltonian function H(q,p ) 
that lead to the collapse of the adiabatic evolution. The global bifurcations correspond to the onset 
of tunnelling in QAA and lead to the failure of the algorithm. In contrast, the local bifurcations 
while still corresponding to exponentially small minimum gap only lead to the decrease of the 
probability of success by a factor of 2. Since in a given problem function H(q,p) is a low de- 
gree polynomial in its arguments there are only a few local bifurcations possible. However, the 
phenomenon of local bifurcations may become important for more difficult random optimization 
problems. Assuming the number of such bifurcations M is large the probability of success is 
reduced by a factor of 2~ M . For M that scales up with n that would lead to the failure of the 
algorithm. 

The method described in the paper, suggests an interesting extension to the case of the K-SAT 
HWP in the limit of large K. In the large-spin limit, the effective potential can be described as a 
white noise with certain intensity and the well-known methods (optimal fluctuation approach) are 
applicable. This leads to the complete characterization of the minimal gap statistics. In particular, 
there always exists a range of parameters where the gap and therefore the run-time of QAA are 
polynomial in n in the large-n limit. 
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APPENDIX A: SPIN OPERATOR REPRESENTATION OF MATRIX A 


The random real symmetric 8x8 random matrix A introduced in (17), describes the ’’tran- 
sitions” between each of the 2 3 = 8 states for each clause involving 3 bits [1]. This matrix has 
(8 x 8 - 8) /2 = 28 independent matrix elements and can be presented in the form 

A = A (1) +A (2) +A (3) , (Al) 

where A^\ A^ 2 \ A^ correspond to the transitions involving one, two and three bits, respectively. 
For each realization, we have 

A^ = a a a x - b ss > (l 4- s a z p) (l + s a*) , (A2) 

A (2) = a a 0 (erf of + a~a^) ^ 6 S (l + s cf) 

S=± 1 

+ a a0 (<rf trj + <T~a+) ^ (1 + s a i) > 

A (3) = B (at at <rf + erf a f erf) + C (at at 'of + erf a f <rf) 

+D (erf at at + af <xf erf ) +E (erf erf cr+ + erf erf erf ) , 

where 

a<y •, bs.s' > a a /3 1 a a0 , b s , B a0 ^ , Coin'll -^a/3'f i E'apy (A 3) 

are the real coefficients. The indices (a,/?, 7 ) € (1, 2,3} label the bits, erf, a x and er* are the 
Pauli sigma-matrices of raising/lowering, ^-projection and ^-projection, respectively and s = ±1 
.is a spin projection variable. Note that the operator | (1 + s erf) is a projector onto the spin state s 
for the bit k. Clearly, the number of independent parameters in (A2) is3x4 + 6x2 + 4 = 28, 
where the three terms of the sum correspond to A^, A® and A^ 3 f respectively. Note that this 
number of parameters equals the number of independent matrix elements of A estimated above. 



24 


In the matrix form, the representation (A2) yields 

0 Q3&-H- (32^++ 0,22 b+ Olb++ <^13^+ Ol2 b+ B 

0 ( 323 ^ 4 - 0,2 bj ^ 13 ^+ ®1^+— C 0,12b— 

0 036+- ai2&+ D aib- + aisb- 

0 E 0 , 12 b- aub- a\b 

0 <236 — p 02b — ). 023b— 

0 023b- 02 b — 

0 a 3 b 

0 




where the vector £ of 8 basis states is 

£ — [+ + +; +H — ; H b; +H — ; 


+ +i 


+ 


+; . 


(A4) 


(A5) 


and the lower left portion of the symmetric matrix A c is obtained by reflection with respect to the 
diagonal. 

The driver He is obtained by summation over all clauses. In doing this summation, we take 
into account that now the bit indices {i,j, k} G {1,2, ...n} run through all n bits, whereas the 
indices (a, (3, 7) G {1, 2, 3} characterizing the realization of Ac still run through the 3 bits (since 
the same realization of Aq is applied to all triples of bits). The driver He is given by 

H e = 4 !) + + Hf, (A 6) 

where Hg\ Hg \ Hg' correspond to the transitions involving one, two and three bits, analogous 
to (Al). 

H [ e = bss> ( 1 + S(jZ j) (l + 5 ^)> (A7) 

a y i f j,k s,s'—±l 

h e ] = J2 a *0 (°t +a i a i) \ bs ^ + s a *) 

a,j3,i,j,k s=± 1 

+o a(3 (at <J- + of <rt) ]■ b * (1 + 5 a k) » 

S=± 1 

h e ] = a t °k + *f a J °k) + C {°i G k + of oj ) 

hi ft 

+ D (erf erf <x+ + of aj erf) +E (a+ o~ c£ + of cr+ af ) . 

One should note that the second term on the r.h.s. for gives a contribution, which is diag- 
onal in S z representation and therefore leads to the effective ”re-definition” of the cost function. 
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Following the logic of [1], we disregard such terms. Also, the commutation relations between the 
total spin components give contributions ~ l/l to the effective potential and can be neglected in 
the large-spin limit. Taking this into account, we obtain from (A 7 ) in the large-spin limit 

/ 7T. \ 3 

h e = ( 7 i + 72 n\ + 73 + 74 n x n z + 7 5 n x n 2 z + 7 6 n 2 x n z ) , 

where n a = S a /l is a dimensionless spin projection on a-axis and the coefficients {7/.} are given 
by 



In particular, the deterministic driver considered in [ 1 ] corresponds to 6_ + = 6 + _ = {&,.} = 

= B = C = D = E = 0,a a = l and b ++ = —b = —2. It follows from (A8) that in 

this case, the only non-zero coefficient in (A8) is 74 = - 8 . This corresponds to H E = — 4 nS x S z , 
which is equivalent to H E = —2 n (S X S Z + S Z S X ) in the large-spin limit according to the above 
discussion. 


APPENDIX B: BIFURCATION POINT ANALYSIS 


Taking into account only the 74 term in H E and expanding up to the 4 th order, we obtain the 
conditions U' — U" = U m = 0 for the A 3 bifurcation point {r c , j c , x} in the form [ 24 ] (cf. ( 52 )) 

3 


2x -x 3 + 7 c r c 1 - -x 2 
\ 2 , 


T c (A + 2Aa; + 3A^ 2 ) = — (1 — r c ) 

r c (2 A + 6At) = - (1 - r c ) [2 - 3 x 2 + 7 c t c (- 3 a:)] , 
6r c A = (1 - r c ) ( 6x -(- 3t c t c ) , 


(Bl) 
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which have to be solved for j c , r c and x for the given {/?*}. Solving for x, we obtain condition 
(53) in the text. 
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